Shifts of Random Energy Levels by a Local Perturbation 
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We consider the effect of a local perturbation on the energy levels of a system described by 
random matrix theory. An analytic expression for the joint distribution function of initial and final 
energy levels is obtained. In the case of unitary ensemble we also find the two-point correlation 
function of initial and final densities of states. 
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The random matrix theory of energy levels in 

complex systems was developed in the fifties for the de- 
scription of the absorption spectra of large nuclei. In this 
approach one gives up any attempt to study the position 
of each particular resonance, but instead concentrates on 
the characteristics averaged over large number of reso- 
nances. The positions of the resonances are identified 
with the eigenvalues of some matrix, and the averaging 
is performed over the elements of this matrix. 

The statistical properties of the eigenvalues Si of a her- 
mitian random matrix H of size N x N are completely 
described by their joint distribution function P{{ei}), 
which can be written in the following simplified form 

P{{e.})^Y[{e^-e,f. (1) 

i>j 

Here the energies Si are ordered: Si < the expo- 

nent /3 = 1,2,4 for ensembles of orthogonal, unitary and 
symplectic matrices, respectively. 

Expression (|^) vanishes when any two energy levels ap- 
proach each other; this effect is commonly referred to as 
level repulsion. This repulsion can also be illustrated by 
rewriting P{{ei}) as 

P({£0) = e-'^^, £;^-^ln(e,-e,). (2) 

Thus P{{ei}) can be interpreted as Gibbs distribution 
of a gas of classical particles at points Si with loga- 
rithmic repulsion between them. In physically interest- 
ing situations there must be a finite average distance A 
between the particles (energy levels). This is usually 
achieved by either introducing a parabolic confining po- 
tential SE = a or by confining the particles to a 
circle [||. 

More recently the random matrix theory was applied to 
a number of physically different systems, such as metallic 
grains j^,^ and microwave cavities 0. In these systems 
one can easily modify the matrix H, e.g., by applying 
magnetic field to the grain or by deforming the microwave 
cavity. Such modification can be described by a pertur- 
bation V, and one is usually interested in the correlations 



of the energy levels of the old and new systems described 
by matrices H and H + V. 

Dyson Q suggested to describe such correlations in 
terms of the viscous Brownian motion of the infinitely 
heavy particles (|^) with the same logarithmic interac- 
tions between them. The external perturbation in this 
approach plays the role of the fictitious time in which 
the Brownian motion occurs. What remains then is to 
find the distribution of the positions of all the particles 
after some time t, provided initial distribution ([^). Time 
t can be related to the characteristic value of the poten- 
tial V, so that the parametric correlations are universal 
functions of only one parameter. 

However the Brownian-motion approach is not always 
applicable. Consider a perturbation of the general form 

N 

V = NJ2v^\^){^l (3) 

i=l 

where |i) form a complete set of states, and the matrix 
dimension N is included for proper definition of limit 
N oo. For the Brownian-motion model to be applica- 
ble H,^, the condition Vi A must hold, where A is the 
mean level spacing. Since ^ 1, the sum J^i ''jf which 
has the meaning of the fictitious time t in the Brownian- 
motion picture, can still be arbitrarily large. 

In a number of interesting physical situations one deals 
with a local perturbation described by Eq. (^ with 
Vi = v5ii, where v is not necessarily small. An exam- 
ple of such perturbation is a short-range impurity in a 
metallic grain. When such an impurity is added to the 
system, its levels Si shift to new positions A^. The new 
many-particle ground state |^) has a rather small over- 
lap with the old one, 1$); this phenomenon ||l^ is called 
orthogonality catastrophe. The overlap can be expressed 
in terms of the old and new energy levels: 

K*i*>i^^n if i^;::!!;:::! - w 

where M is the number of electrons in the system [ pT| . 
Therefore to treat the orthogonality catastrophe in a 



1 



metallic grain one needs the knowledge of the joint distri- 
bution function P{{ei}, {Xi}) of both old and new energy 
levels. 

The determination of this joint distribution function 
is the main subject of this paper. The Brownian-motion 
model is not applicable in this case, but a closed analytic 
expression for P({ei}, {Xi}) can be found directly. For 
the orthogonal, (3=1, and unitary, /3 = 2, ensembles we 
will show that 

P{{e^},{X^}) QC 



X exp 



(5) 



Energy levels in Eq. (o) are constrained by the condition 



Ei < Xi < Si+i, if w > 0, 
< Xi < Ei, if V < 0. 



(6) 



Equation (||) is the central result of this paper. 

To derive Eq. we need to relate the eigenvalues Xj 
of the perturbed matrix H + vN\v){v\ to the unperturbed 
eigenvalues Si and eigenfunctions here \v) is an arbi- 
trary vector. This is easily accomplished: 



^-^ Xj — Si vN 



(7) 



Equation (0) enables one to find P{{ei}, {Xi}) given the 
joint distribution function of unperturbed eigenvalues 
and eigenfunctions. Since the distributions of eigenval- 
ues and eigenvectors in the random matrix theory are 
uncorrelated, we have 



P{{e.},{^^})^Pi{e^})p{{A.}) 



det 



dA,, 



dXi 



(8) 



Here P{{si}) is the distribution function (|l]) of the imper- 
turbed energy levels, and p({Ak}) is the eigenvector dis- 
tribution function, which at A'^ — > cx) is given by Porter- 
Thomas ilSf formula: 



N 



'2nNA,)^~P/- 



■ exp 



NA, 



(9) 



where /3 = 1 or 2. Finally, the last factor in Eq. ^ is 
the Jacobian of the transformation from the eigenvector 
variables Ai to the new energies Xj . 

To find the distribution function (H) we do not need 
to solve the equation with respect to Xj. In order 
to find the Jacobian in Eq. (^) one only has to solve (j^) 
with respect to Ai. The latter is a much simpler problem 
since equation (Q) is linear in Ai, and the solution can be 
expressed in terms of Cauchy determinants. This readily 
yields 



1 nf.jA, 



e^) 



vN Urn 



(10) 



By definition all Ai are positive, see Eq. (0). This im- 
mediately gives constraint (H). It follows from Eq. (|lo| ) 
that 



dX, 



A 



X,-e 



As a result the Jacobian in Eq, 
determinant, and we obtain 



det 



dA^ 
dXi 



1 



(11) 

is reduced to a Cauchy 



(12) 



In order to express the Porter-Thomas distribution 
function (|p) in terms of the energies Ei and Xj we need 
to evaluate the sum Ai . To this end we sum up both 
sides of identity (ITTI ) over i, and using Eq. d), find 



J - 1, 



We therefore conclude that the sum of Ai is a linear func- 
tion of all Ai. The constant can be determined by noticing 
that according to Eq. (|o|) at — > Si we have Ai ^ 0. 
Thus 



Nv 



(13) 



Finally, we substitute Eqs. (0), (0), and (|2|) into 
Eq. d), and with the help of Eqs. ^ and (|l|) get the 
joint distribution function Eq. Strictly speaking the 
result (H) is valid only in the limit — * oo. This is the 
physically most interesting regime where the properties 
of the system are universal. 

If one is interested in non-universal corrections associ- 
ated with the finite size of the matrix, the Porter-Thomas 
distribution (^) should be replaced by [Q 



p/2-1 



S[l-J:^A.]■ (14) 



As a result, the exponential factor in Eq. (||) is replaced 
by S[l — {Nv)~^ J^ii'^i ^ ^i)]- However, for practical cal- 
culations the distribution function in the from (^ is more 
convenient. 

Our result (^ contains complete information about 
distribution of the old and new energy levels of the sys- 
tem. In applications, such as analysis of experimental 
spectra, one often needs only a small part of this informa- 
tion, which is contained in n-level correlation functions. 
The most important of them is the two-point correlation 
function of the old and new densities of states K2. We 
define this quantity as 



2 



ik 



Sk 



K2{s) has the meaning of the probabiHty to find a new 
level at a distance s from a given old level. The corre- 
lation function K2 does not depend on energy E, pro- 
vided that s is much smaller than the characteristic en- 
ergy scale over which average density of states p{E) — 
i^- ^i)) varies. 
We have been able to obtain a compact analytic ex- 
pression for K2 (s) for the unitary ensemble only, and we 
outline the derivation below. Because K2 does not de- 
pend on the particular shape of p{E), it is convenient 
to get rid of energy dependence of p{E) by adopting the 
circular ensemble of Dyson ||^ , where all TV levels are put 
on the circle of unit radius. Equation (||) for (3 = 2 then 
takes the form 



sm 



sm 



Ai-A, 



(15) 



Here the energies are measured in dimensionless units 
and defined within the interval [— tt, tt] ; we have also omit- 
ted the normalization constant. Mean level spacing in 
such model is given by A — 2tt/N . 

We express function K2 in terms of the functional 
derivative 



K2 - 



A2 5^1 [A, B\ 



I [A, B] 5A5B 



(16a) 



A.B=0 



of the generating functional for u > 0, (negative v are 
considered analogously) 



1=1 dei / d\i / de2 / dX2 ■ 
l[{l + A{ek)){l + B{Xk)) 



dX^P, 



(16b) 



Following the procedure similar to that of Ref . y| , we 
find that the generating functional / can be rewritten as 
determinant of certain matrix / = det where 

Fki = Cde CdX e'^^-'+^k)e-(v-^+a)\ ^ ^^^^^j ^ 

J — IT J e 

(17) 

We then expand the determinant up to the second order 
in small sources A and B. This expansion requires the 
knowledge of the matrix a.t A = B = 0, which in 
the hmit ^ 00 takes the form FfTi^ = 5ki " It''^ ■ Sub- 



stituting the result in Eq. (16a) we obtain after simple 
algebra [O 



dr'e " 



A sm nr 



d_ 

dr 



rA sm nr 



(18) 



where r = s/A is the energy in the units of the level 
spacing. The result (18) is valid for positive v; for w < 
one should substitute r — r in Eq. (p^). 

Let us now discuss asymptotic behavior of the two- 
point correlation function (|l8[). In the limit of vanish- 
ing perturbation w — > 0, we immediately obtain K2 = 
6{r) -|- 1 — (sin7rr/7rr)^ which is the well-known result 
for the two-point correlation function of the unitary en- 
semble H:^. In the limit r ;2> l,w/A, the integration in 
Eq. (|l8| ) can be easily performed and we find 



K2 = l- 



sin^ (irr — S) 



(19) 



where 5 = arctan (ttw/A) is the phase shift of the scatter- 
ing off of the impurity. The sequence of periodic maxima 
of the correlation function ([l9| ) is a signature of the level 
repulsion, and the average shift d of the new levels with 
respect to the old ones is consistent with the Friedel sum 
rule (A, - e,)/A = S/n. 

In conclusion, we studied the statistics of shifts of 
eigenvalues of a random Hamiltonian by a local pertur- 
bation of arbitrary strength. Despite the fact that the 
conventional Brownian motion model is not applicable, 
we have found the whole joint distribution function of old 
and new levels, Eq. (||) , and the correlator of the old and 
new densities of states (|l8|). Using Eq. (^), these results 
can be applied to the study of orthogonality catastrophe 
in small metallic grains. 
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